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We have studied the weakly non-linear quantum transport properties of a two- 
dimensional quantum wire which can be solved exactly. The non-linear transport 
coefficients have been calculated and interesting physical properties revealed. In 
particular we found that as the incoming electron energy approaches a resonant 
point given by energy E = E r , where the transport is characterized by a com- 
plete reflection, the second order non-linear conductance changes its sign. This 
has interesting implications to the current-voltage characteristics. We have also 
investigated the establishment of the gauge invariance condition. We found that 
for systems with a finite scattering region, correction terms to the theoretical 
formalism are needed to preserve the gauge invariance. These corrections were 
derived analytically for this model. 
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1 Introduction 



Non-linear quantum transport in mesoscopic systems has been a very active re- 
search field in recent years 0, |3|, J|. Taboryski et. aZ.[Q have reported obser- 
vations of non-linear and asymmetric conductance oscillations of quantum point 
contacts at a small bias voltage. They found that the non-Ohmic and asymmet- 
ric behavior causes a rectified dc signal as the response to an applied ac current. 
On the theoretical side, Wingreen et. a/.@ have presented a general formula- 
tion to deal with the situation of a non-linear and time- dependent current going 
through a small interacting region where electron energies can be changed by 
time-dependent voltages. At the same time, Biittiker and his co-workers ]|, [l], H 
have advanced a current conserving theory for the frequency dependent trans- 
port. Recently, this current conserving formalism has been applied to a two 
dimensional mesoscopic conductor^]. This theory can also be applied to discuss 
the non-linear behavior of mesoscopic samples and the theory is gauge invariant. 
It has been recognized |§ that in non-linear coherent quantum transport, it is 
essential to consider the internal self-consistent potential in order to satisfy the 
gauge invariant condition. This condition demands that all physical properties 
predicted by a theory can not change if there is a global voltage shift. Obviously 
this is a fundamental requirement. 

Recently, Christen and Biittiker || have investigated the rectification coeffi- 
cient of a quantum point contact and the non-linear current- voltage characteristic 
of a resonant level in a double barrier structure using the theory of gauge invari- 
ant non-linear conductance. Another important application of this theory is to 
investigate two-dimensional (2D) mesoscopic and ballistic quantum devices which 
can now be routinely fabricated in many laboratories. Unfortunately due to a 
particular technical difficulty, namely the evaluation of a quantity called sensi- 
tivity (see below), so far little is known for the non-linear conductance in 2D. 
Clearly an understanding of 2D situation is very much needed in order to gain 
further intuition to the coherent transport and to predict the non-linear char- 
acteristics of the variety of 2D nanostructures. The purpose of this paper is to 
investigate the gauge invariant non-linear transport in a specific two dimensional 
system which can be solved exactly. Hence we were able to obtain various rel- 
evant physical quantities. Although a general study for an arbitrary 2D system 
seems difficult, our perspective is that an exactly solved model is valuable since 
it clearly and unambiguously reveals the physical properties of the non-linear 
transport coefficients. 
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To be specific, we have considered a very simple 2D model which is a quasi-lD 
ballistic conductor with a ^-potential confined inside, as shown in figure (la). 
Because quantum scattering in this system leads to mode mixing which is the 
basic feature of a two-dimensional system, it provides answers to our 2D prob- 
lem. In a previous work(13| we have used this model to study the electric current 
conservation of the AC-transport formalism at the linear conductance level, and 
calculated the important physical quantities such as the global and local partial 
density of states. In the following we shall extend our calculation to explicitly 
calculate the second order non-linear conductance G m and Cm- Due to the 
gauge invariant condition (see below), we should have G m + Gu 2 = 0. It turns 
out that for systems with a finite scattering volume, as those of any numerical 
calculations, we found that a correction term must be added to satisfy this con- 
dition. For this system there is a resonant state with energy E r characterized by 
a complete reflection, i.e. the reflection coefficient R — 1. Our results showed 
that the second order non-linear conductance G U1 changes sign near the resonant 
point E r . This leads to interesting current- voltage characteristics of this system. 

The paper is organized as the following. In the next section we shall briefly 
review the gauge invariant theory for non-linear transport set out by Biittiker[|l|. 
In section III we will present the solution of the 2D scattering problem. Some of 
the technical details of section III have been put into the Appendix. Our results 
are presented in section IV. The last section serves as a brief summary. 



2 Gauge invariant formalism 

To be complete, we shall first briefly review the gauge invariant formalism of 
Christen and Biittiker[|| and set out our calculation procedure for the 2D system. 
For a multi-probe mesoscopic system, the current through probe a is given by]]], || 

J « = T £ / dEf ( E ~ Ef ~ eV P) A «p( E ' TO). (!) 

p 

where f(E) is the Fermi distribution function, and 

A a p(E, {V 7 }) = Tr[l a 5 a(3 - slp(E, {V 7 })s a p(E, {V 7 })] (2) 

are the screened (negative) transmission functions. For the weakly non-linear 
transport, Eq.([l]) can be expanded with respect to the voltages Vp, 

i a = G a pVp + G a pjVpVry + (3) 

p p<r 
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where 

2e 



G aP = -^J dE{-d E f)A aP (4) 
is the linear conductance and 
e 2 f 

G a f3-y = -r J dE(-d E f)(d Vl A al3 + d VfJ A ai + ed E A aP 5^) (5) 

is the second order non-linear conductance. In Eqs.((|) and fl5|), the A a p are 
evaluated at {V^} = 0. The requirements that the current is conserved and be 
independent of a global voltage shift (gauge invariance) yield [|l], j|, |§ 

G a f3 = X! G a p = 

and 

a (3 7 

From this equaiton and Eq.fl5|), the gauge invariance condition for A a p is 

ed E A a p + J2d Vj A af i = . (6) 

7 

The derivative dy^A^ can be expressed in terms of functional derivative of A a p 
with respect the electric potential U and the characteristic potential u 1 which 
satisfies J2^ u -y — 1 

d Vl A aP = J dh^u^r) . (7) 

To gain further insight on Eq.(|6|), let's consider a two probe system. Eq.([3|) 
can be written as 

h = G n V! + G l2 V 2 + G m V? + 2G 112 V 1 V 2 + G 122 Vi 

Obviously, G% 2 = —G\\ due to the conservation of electric current. From Eqs.@, 
(||), and (0) and using the fact that U\ + u 2 = 1, we have G m = —Gn 2 = G\ 22 . 
Therefore the current depends only on the voltage differences which is the direct 
consequence of the gauge invariant condition Eq.@. We obtain, 

h = G 11 (V 1 -V 2 ) + G U i(V 1 -V 2 ) 2 . (8) 

To calculate the transmission functions A a p and their functional derivatives, 
we need the characteristic potential u 1 . This, in turn, needs the solution of 
the Poisson equation with a nonlocal screening term[|l], ||. To actually carry 
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out this procedure is very complicated. However if we can use the Thomas- 
Fermi approximation, which is more appropriate for metallic conductors, the 
characteristic potential is simplified and is found to be related to the local partial 
density of states. Within Thomas-Fermi screening, we obtain 

_ dn(r, 7 ) dn(r) 
%l } ~ dE 1 dE ' W 

where the partial local density of states dn(r, r y)/dE is called the injectivity and 
is given by|| 

jraQvr) _ l^ 7 nl 2 (10) 

dE ^ hv jn ' 1 ' 

where v in is the channel velocity and \l/ 7n is a scattering state. Finally the term 
dn{r)/dE = J2 a dn(r,a)/dE is the total local density of states. Substituting 
Eqs.(|]) and (§) into Eq.(|7|), we obtain 



where 

^ = s lx8TFT7Z\ + S <*PT777^\ ( 12 ) 



" a nU{r) +Sal3 5U{r) 



is called sensztm%[|10| . We are aware of two ways of calculating the sensitivity^]. 
The first is to evaluate 5s a/ 3/SU directly by introducing a (^-function of infinites- 
imal strength SU inside the scattering region. Alternatively, one can calculate 
it using the retarded Green's function. For a 2D system, in general the Green's 
function can not be obtained explicitly, hence we shall use the first method by 
directly compute the sensitivity. After obtaining the sensitivity, we can then 
compute dv 1 A a p from Eq. (|ll|), and obtain G Q/ g 7 from Eq. (|5|). 



3 Model and Analysis 

As mentioned in the introduction, figure (la) shows the system where a 5- 
potential is confined inside a quasi-lD wire with width a. We assume, for sim- 
plicity of the calculation, that the boundaries of the ballistic conductor are hard 
walls, i.e. the potential V = oo at the walls. Inside the conductor, the potential 
is zero everywhere except that a 5 function potential V(x,y) = •y6(x)5(y — yo) is 
placed at position r = (0,y ). The scattering region x\ < x < x 2 is assumed to 
be symmetric with x 2 = —x\ = L/2. From now on we set H = 1 and m = 1/2 to 
fix our units. 
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The transmission and reflection amplitudes have been calculated using a mode 



matching method[12, O. When the incident electron is in the first subband, in 



an earlier work we have explicitly obtained these amplitudes [13]. The evaluation 
was straightforward but quite tedious, we refer interested reader to Ref. |H| for 
details of this algebra. Here we only quote the results: for reflection the amplitude 
is 

b n = (13) 

and for transmission it is 

c„ = S nl + b n . (14) 

Here a = 1 + i £ n T nn /(2k n ); T nm = 7X%o)Xm(2/o), and Xn{y) is the wave 
function of the n-th subband along y-direction. k n = E — (nir/a) 2 is the longitu- 
dinal momentum for the n-th mode; % = \/— 1. Note that for electron traveling 
in the first subband, k n with n > 1 is purely imaginary. For our coordinate 
system the scattering matrix elements s a p are given by sn = biexp(ikiL) and 
S12 = Ci exp(ikiL). 

As mentioned in the last section, to calculate the non-linear conductance of 
our 2D sample, it is necessary to find the sensitivity r] a p. Hence according to 
Eq. (|i~2| ) we must evaluate Ss a p/SU(xi,yi) where the pair (xi,yi) is an arbitrary 
location in the scattering volume. For a general 2D sample a direct calculation 
of this functional derivative is very difficult if not impossible. Fortunately for our 
model this can actually be done exactly. As a first step we shall introduce an 
additional (^-potential of infinitesimal strength SU at position (xi,yi) inside the 
scattering volume. Thus our system becomes to that shown in figure (lb). Then 
we shall solve the scattering matrix formally as a functional of SU. Obviously 
being able to carry out this step is crucial. Finally the functional derivative is 



performed. To proceed we again use the mode matching method [12, 11, 13]. We 
will assume x\ < in the following calculation. The calculation for x\ > can 
be done in a similar fashion. The electron wave functions are written as follows. 
For region I (see figure, (lb)): 



= J2Xn(y)(a n e iknX + b n e 



•iknX \ 



where a n is the incoming wave amplitude and taken as an input parameter; b n is 
the reflection amplitude. Similarly for region II 

*n = Y,Xn(y)(e n e ik ^ + f n e- ik ^) , 
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and for region III 

^iii = Y.Xn(y)c n e lknX , 

n 

where c n is the transmission amplitude. We shall match the wavefunctions and 
their x-derivatives at the positions x — x\ and x — 0. We obtain, at x — x\ 

a n e iknXl + b n e~ ik " Xl = e n e iknXl + f n e- iknXl , (15) 

and 

irvnydnfi J n& ) lrun\ u "n t ' "nt ) 

_ f (p P ik mxi , r -ik m xi\ 

/ j L nm\ c -m { -- ~ J J 

m 

(16) 

where f nm = 5Uxt(yi)Xm(yi)- 
At x = 0, the matching gives 

and 

ik n C n ik n {(i n fn) ^ ^nmC-m 

m 

To simplify the notation, from now on (xi, yi) is replaced by (x, y). From the last 
two equations, we solve for e n and /„ 

2?/c n 6 n ^ ] Pnrr&m i ( 1 7) 

and 

2ik n f n = ^ , r ^ nrrfim, ; (^) 
m 

where P ram = T nm — 2ik n 5 nm . Eliminating 6 ra from Eqs.(15) and (|l6|), we obtain 



9iJc p p iknX — Oih n p iknX J-Vr P ifc ™ x _|_ f p-'^ 1 ! 

m 

Taking the limit 5U — > 0, we have 

2 ^-T77 = E F n n (fl w e tt ^ + b m e~ ik ™ x ) , (19) 
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where Y nm = T nm /5U. To arrive at the above result we have used the fact that 
as 5U —>■ 0, i.e. when the extra 5 function vanishes, we must have e m = a m and 
f m = b m . From Eq.(0), 



n 5f/" ^ nm SU- 
From Eqs. ([l9|) and (p0|) , we arrive at 

V f (a e ikmX + b e~ 



(20) 



_ V P —e" 

/ j 1 nm. ~ , ' 



6U 



(21) 



where \? = ^fj for x < 0. From Eq.(21), we have 



5cj, 
5U 



^ = -£(p-%e-^x;* • 



The matrix P 1 has been obtained in Ref. [13|, and we quote 

From this equation and Eq.(|l4|) we see that for I = 1, i.e., the first subband, 
(P ^in = ic n /(2ki) provided that \n is real which is true in our case. This 
yields 

8c\ 



Similarly, from Eqs. QI^), (0), @, @, and (0), we obtain 



(22) 



6K 
5U 



" e n 2ik n x i 

p 2ik n x x„ 1 r„ 



2zAv, 



sin(fc n x 



5c n 



When n — 1, Sbi/SU becomes, 

<H>1 1 



<H7 



2zA;i 
1 

2i£^ 



(Xie ifcia ' + ^6nXne" 

^ 2 



(23) 
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Because the scattering matrix elements sn ~ b\ and S12 ~ c\ as mentioned 
above, with the functional derivatives Eqs. (^) and ( |22|) we can evaluate 5s a/ 3/SU 
trivially thus obtaining the sensitivity 7] a p of Eq. (|T2|). Then using the prescrip- 
tion discussed at the end of Section 2, we obtain all the weakly non-linear con- 
ductances and other quantities of interest. Our results will be presented in the 
next section. 

To end this section of the theoretical analysis, we mention that to check the 
result of functional derivatives, in the Appendix we shall explicitly calculate a 
quantity called emissivity^ using these functional derivatives. In the absence of 
a magnetic field, it is known|§ that emissivity equals to the injectivity defined 
in Eq.(|l(]) which we can compute using the wavef unctions. Indeed we confirm 
in the Appendix that these two equal thus providing a necessary check to the 
calculations presented here. 



4 Results 

To obtain numerical results from our analytical formula, for the system of figure 
(la) we consider incident electron coming from probe 1 and set a — L — 1, 
yo = 0.3, and 7 = — 1. Although we have restricted the incoming electron energy 
to the first subband, quantum scattering at the (^-function potential leads to mode 
mixing. Thus in our numerical calculations we have included 50 modes in the 
scattering volume. We have checked that this is enough to obtain good numerical 
convergence. 

As a first result we plot the sensitivity rju (r, E) as a function of the electron 
incident energy E at several positions r. This is shown in figure (2). As discussed 
in Section 2, r] a p appears naturally in the theoretical formalism, and it essentially 
describes the local electric current response of the scattering problem when there 
is a small local potential change. It is related to the real part of the diagonal ele- 



ments of the Green's function JT0|. figure (2) not only shows interesting behavior 



of this quantity, but also gives vivid intuition about the local current response. 



As shown in our earlier work||13|| and mentioned above, the quantum wire studied 
here has a resonant state at energy E = E r = 36.65 where we have a complete 
reflection (reflection coefficient R — 1). Not surprisingly, the sensitivity has a 
large peak at this resonance energy because the system response is most sensitive 
to potential perturbations at resonance. On the other hand, this peak is larger 
when we are closer to the (^-function scatterer located at x = 0: this indicates 
that the local perturbation has larger effects when it is closer to the scattering 
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center. Although figure (2) shows r/u at positions to the left of the scatterer, 
we have checked that its behavior is exactly the same for positions x > as our 
system is symmetric. 

Adding up all the local responses according to Eq. (|11]), we can explicitly 
examine the gauge invariant condition Eq.(|]). Using Eqs.(0) and (0) and the fact 
that u\ + U2 = 1, Eq.(||) reduces to 

It is straightforward to evaluate the left hand side of this equation. Using the 
functional derivatives obtained in the last section, as well as the energy derivatives 
of Eqs. fll3"l) and flip), we found that the left hand side of the above equation is 
nonzero, and is given by 

correction = ^&(s n ) + i?e(E ^r~(e iknL ~ 2 )) • (24) 

Thus in order to have precise gauge invariance, this correction must be included. 
From this result, we notice that the first correction term is only significant near 
the first subband threshold where k\ ~ and is negligible for larger incoming 
electron energies. For the second correction term, let's examine its behavior near 
the n-th subband with n > 1 . From Eq. fll3"l) we see that as the incoming electron 
momentum k — ► k n , h\ — * k n and b n is finite. Therefore, the second correction 
term remains finite when electron energy approaches the n-th subband (n > 1). 
This is different from the AC-transport where the correction diverges ||13|| near 
the n-th subband with n > 1. We emphasize that the correction term is due 
to the fact that we are considering a finite scattering volume. As the scattering 
volume or the incident energy become larger, the effect of these correction terms 
diminishes. This can be seen clearly due to the factor k\ in the denominator, 
and the exponentially decaying factor exp(i/u n Z/) as k n is purely imaginary for all 
n > 1. 

Now we present numerical evaluations of the analytical formula derived in 
the last section for the second order non- linear conductance G a ^. In figure 
(3) for a comparison the solid line shows the linear conductance G21 which is 
proportional to the transmission coefficient by the Landauer formula. The dotted 
line shows the second order non-linear conductance G\\\. In the presence of an 
attractive 5-function scatterer, there exists a quasi-bound state at E = E r = 
36.65. As a result, we observe the complete reflection 0, |l3j indicated by G2i=0. 
As expected, the non-linear coefficient G in also vanishes at E r . Furthermore 
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Gin changes its sign as the incoming electron energy passes through the resonant 
energy. This has important implications on the current- volt age characteristics if 
we recall the I-V relation Eq. (^). The I-V curves of this quantum wire system 
is shown in figure (4) for several different electron energies. We can clearly see 
that when E is smaller than E r , e.g. for E = 10.91 and 36.22, since both G 2 i 
and Cm are positive, the current I± increases with potential difference of the 
two probes AV = V\ — V%. However when E = 37.11 which is just above E r , I\ 
is a decreasing function of AV due to the negative non-linear conductance G m . 
For this case the quantum wire has a negative differential resistance. Finally at 
a even larger energy E = 39.03, while G m < 0, the value of G 2 i is large enough 
such that the linear contribution dominates at small AV and the non-linear term 
is larger at larger AV. This behavior is shown in the dash-dotted line in figure 
(4). 

To numerically demonstrate the gauge invariance, we have also computed the 
non-linear conductance Gn 2 according to Eq. @ and it is shown in figure (5) 
as the dotted line. For comparison we have re-plotted the non-linear conduc- 
tance Gm (solid line). Note that near the resonant region, these two non-linear 
conductances are very close to each other. This is a surprising result, because if 
the gauge invariant condition is precisely satisfied we should have G m = — Gn 2 . 
Thus without the correction term discussed above, the gauge invariance could 
not even be satisfied up to a sign ! Hence it is important, in any numerical calcu- 
lations, to have a large scattering volumef^]. Indeed, when we add the correction 
term of Eq.(^4j) to G112, (dashed line in figure (5)), we obtain the expected perfect 
agreement to the gauge invariance. 

5 Summary 

To summarize, we have solved exactly the weakly non-linear transport character- 
istics of a two-dimensional quantum wire model. To the best of our knowledge 
this is the first exact solution for a truly two-dimensional ballistic model. The 
second order non-linear conductances are derived analytically. We found that as 
the incoming electron energy crosses the resonant point, the non-linear conduc- 
tance changes its sign. This leads to interesting current-voltage behavior when 
the incoming electron energy changes. We have also examined the gauge invariant 
condition which is obtained by the global voltage shift. We found that for sys- 
tems with a finite scattering volume, correction terms are needed to preserve the 
gauge invariant condition. We have derived these corrections analytically for our 
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model. The correction term consists of two parts. The first part dominates when 
the incident energy E is near the first subband threshold. On the other hand 
the second part is given by the amplitudes of the non-propagating modes and 
is significant near the resonant point. Finally, our exact calculation reveals the 
interesting behavior of the sensitivity which describes the local electric current 
response to a potential perturbation. 
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Appendix 

To check our result of the functional derivatives, i.e. Eqs. (|22"D and (p3|), in this 
Appendix we compute the emissivity defined as|T[ 

dn(a,r) L^Tt \ t SSa P _ Ss °P 1 

dE ~ 4vrzY [SaP 5U{v) 5U{v) SM ' 

It has been shown that in the absence of a magnetic field the emissivity is equal 
to the injectivity defined in Eq . (|T0|) . We shall explicitly perform the functional 
derivatives to confirm this fact, hence provide the necessary check to our algebra. 
Using Eqs. (|22|) and (^), we have 



S„ hS 19 = Ci h Oi 

11 5U 12 5U 1 5U l 5U 



2ik 1 
1 



(h* lX ie iklx + Xie~ lklX + (1 + 26*) ^ KXneT*™)* (25) 



n=2 



where the relation c\ — 1 + b\ has been used. Before we proceed further, let 
us derive a useful relation from the unitary condition of the scattering matrix, 



12 



namely 

1 + 26? = -T 1 = — . 26 

bi a* 

The first equality comes from the unitary condition c\bi + cib\ = or b\ + (1 + 



26*) bi = 0; and the second equality is from Eq.(13). Since the incoming electron 



is in the first subband, we have fc* = — k n for n > 1. Hence for n > 1, 

K a 

Substituting Eqs.(^) and (|27|) into Eq. fl25|) , we obtain 

which is equivalent to Eq. fllCf ). Notice that the imaginary part of left hand side 
of Eq.(^) is proportional to the emissivity. Its real part gives the sensitivity 
Vn + Vi2- From the unitary condition we have rju + r] 12 = which agrees with 
Eq.©. 
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Figure Captions 



Figure 1. Schematic plot of the quantum wire system, (a). The quantum wire system 
we have studied: a 5 potential j5(f — fo) is confined inside a quasi-lD 
quantum wire, with r~o = (0, yo)- The wire width is a. The scattering 
region is between x\ and X2, where X2 = —X\ = L/2. In our calculations, 
the parameters are set to L = a = 1, y = 0.3, and 7 = —1.0. (b). To 
compute the functional derivatives of the scattering matrix with respect to 
a local potential change, we add another 5 function potential at the position 
(xi,yi). In this case the system is divided into three regions by the dotted 
lines for the boundary matching solution of the Schrodinger equation. 

Figure 2. The sensitivity rju (r, E) as a function of energy E at three different positions 
x = —L/2, — L/4, with the same y = 0.5. For different y the curve r/n as 
a function of E will be multiplied by a constant. Other system parameters 
are the same as those of figure (1). Here the unit of energy is h 2 /(2ma 2 ). 

Figure 3. The conductances G21 and Gm as functions of energy E. Solid line: G21, 
dotted line: G m . Other system parameters are the same as those of figure 
(1). Here the unit of energy is h 2 /(2ma 2 ). 

Figure 4. The current-voltage characteristics as calculated from Eq. (§) at several 
different electron energies E = 10.91,36.22,37.11,39.03. AV = V l - V 2 . 
Other system parameters are the same as those of figure (1). 

Figure 5. A numerical check of the gauge invariant condition. Solid line is G U1 , 
dotted line is G112, dashed line is the G112 + correction where correction is 
given by Eq. (p4|) . Clearly G 112 + correction = -G m . Here the unit of 
energy is h 2 /(2ma 2 ). 
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TempleGraph 2.4, Origin: /morgoth/users/guo@celeborn - Sun Sep 15 19:14:58 1996 (data file was iv_dat) 



Figure (4) 




TempleGraph 2.4, Origin: /morgoth/users/guo@celeborn - Mon Sep 16 09:32:45 1996 (data file was gauge_dat) 



